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ABSTRACT 


This  report  discusses  methods  of  computing  angle  of  attack  by  inference 
using  combinations  of  data  from  presently  available,  on-board  sensors,  thereby 
eliminating  the  need  for  external  vanes  or  probes.  Equations  were  derived  from 
which  computed  angle  of  attack  could  be  extracted.  Those  equations  that  were 
impractical  from  a  mechanization  standpoint  were  eliminated,  leaving  three  can¬ 
didate  methods.  These  three  methods  were  then  analyzed  with  respect  to  errors 
arising  from  mathematical  simplifications  and  errors  due  to  imperfect  sensor  in¬ 
formation.  Two  of  the  candidate  methods  provide  inertial  angle  of  attack,  and 
will  provide  acceptable  accuracy  for  low-performance  aircraft  applications.  The 
third  method  will  provide  a  high-quality,  air-mass-related  angle  of  attack.  A 
mechanization  of  a  high-quality  angle  of  attack  system  using  a  small,  special 
purpose  digital  computer  is  described.  Using  measurements  of  normal  accelera¬ 
tion,  longitudinal  acceleration,  elevator  position,  flap  position,  throttle  pos¬ 
ition,  airspeed,  Mach,  and  dynamic  pressure,  the  system  will  provide  a 
high-quality  angle  of  attack  measurement  applicable  to  any  high-performance  air¬ 
craft  and  competitive  with  current  vane  and  probe  transducers. 
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SECTION  I 
INTRODUCTION 


True  airframe  angle  of  attack  finds  many  applications  as  an  input  to  the 
sophisticated  systems  on-board  the  latest  generation  of  aircraft.  It  has  been 
used  for  weapons  guidance  computation,  autopilot  functions,  speed  command,  stall 
warning,  and  gust  alleviation.  This  multitude  of  uses  requires  an  accurate, 
reliable  measurement  of  angle  of  attack  (c).  Unfortunately,  the  techniques  for 
measuring  a  have  not  progressed  far  oeyond  the  Wright  brothers’  vane  transducer 
used  in  the  first  powered  flights. 

Externally  mounted  angle-of- attack  transducers  have  had  a  rather  poor  his¬ 
tory  of  reliability  and  accuracy.  To  be  effective,  the  device  must  be  mounted 
in  the  free  airstream  to  minimize  fuselage  flow  effects.  This  necessarily  ex¬ 
poses  the  transducer  to  all  the  severe  environmental  conditions  experienced  by 
the  aircraft.  Human  error  has  also  made  a  considerable  contribution  to  a  trans¬ 
ducer  unreliability.  Most  mounting  positions  have  made  the  transducer  a  conven¬ 
ient  handle  or  step  fcr  ground  crews  and  pilots. 

The  inaccuracies  exhibited  by  vanes  and  probes  can  be  attributed  to  the  ex¬ 
treme  range  of  environment  over  which  operation  is  required  and  to  static  errors 
due  to  the  mounting  position.  Sophisticated  approaches  to  solving  these  prob¬ 
lems  have  produced  some  useful  results,  but  at  considerable  expense. 

The  major  drawbacks  of  externally  mounted  transducer^  can  be  overcome  by 
computing  angle  of  attack  from  combinations  of  internally  mounted  sensors.  The 
concept  of  computed  a  Is  not  a  new  one  -  it  had  been  proposed  as  early  as  20 
years  ago.  However,  it  has  never  received  a  thorough  evaluation,  particularly 
in  light  of  the  highly  improved  computing  techniques  now  available. 

This  technical  report  presents  the  results  of  a  studv  of  computed  a.  and  the 
mechanization  thereof.  The  principal  objectives  of  the  study  are  to: 

•  Derive  equations  that  can  be  used  to  compute  angle  of  attack  from 
existing,  internally  mounted  sensors  (Section  II) 


•  Select  those  methods  (equations)  which  appear  to  be  practical  and  determine 
the  errors  that  occur  when  simplifying  assumptions  are  made  (Section  III) 


•  Determine  the  effect  of  sensor  errors  on  the  a  computation  when  using  the 
practical  equations  (Section  IV) 


•  Determine  the  mechanization  of  an  optimum  system  for  computing  angle  of 
attack,  using  one  or  more  of  the  methods  generated  and  refined  under  the 
first  three  objectives  (Section  V).  The  system  should  maintain  0.1-degree 
accuracy  through  =7 5  degrees  of  pitch  and  roll,  =10  degrees  of  side-slip, 
and  +8g  and  -5g  of  acceleration. 
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SECTION  II 

DERIVATION  OF  EQUATIONS 

1.  BASIC  DEFINITIONS 

Before  proceeding  to  the  derivations,  an  explanation  of  axis  systems  would 
be  expedient.  Aircraft  equations  of  motion  are  given  in  terms  of  a  body-fixed 
axis  system  with  its  origin  at  the  aircraft  center  of  mass.  The  x  and  z  axes  in 
the  plane  of  symmetry  are  oriented  such  that  the  Euler  angles  locating  the  body- 
fixed  frame  are  those  angles  indicated  by  the  vertical  gyro,  0  and  0,  shown  in 
Figure  1. 

From  Reference  1.*: 


Ax„„ 

eg 

=  U  +  QW 

-  RV  + 

g 

sin 

0 

(1) 

A*  „ 
eg 

=  X(«)  + 

X(/3)  + 

>T  65 

1 

(2) 

Ay _ 

eg 

=  V  +  RU 

-  PW  - 

g 

cos 

0  sin 

0 

(3) 

Ay„„ 

J  eg 

=  v  + 

6r  r 

(4) 

eg 

=  W  +  PV 

-  QU  - 

g 

cos 

0  cos 

0 

(5) 

Az„„ 

eg 

=  Z(a)  + 

Z. 

6e  e 

+ 

6t 

6T 

(6) 

I  I  -  I  I 

P  -  -p  R  +  Z  QR  -  PQ  =  Lg  ft  +  L  P  +  L  R  +  6 


b  r 
r 


.  I-  -  I 

^  +  ~  I. 


- -  PR  +  -f-7-  (P2  -  R2)  =  ^  C'  +  Mq  Q  +  M6  be  +  6t 


(7) 


(8) 


R 


xz 


♦  i. 


PQ  +  QR  =  N(j8)  +  N  R  +  N  P  +  N 

* 


■V 


(9) 


Stability  derivatives  with  respect  to  speed  are  omitted  because  this  effect 
is  accounted  for  through  use  of  the  varying  dynamic  pressure  implicit  in  the 
dimensional  stability  derivatives.  Details  of  the  aircraft  simulation  are  pre¬ 
sented  in  the  Appendix. 


♦See  list  of  references. 
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Care  must  be  taken  to  distinguish  between  velocities  determined  with 
respect  to  an  inertial  frame  and  air  that  is  moving  with  respect  to  an  inertial 
reference  frame.  The  former  quantity  is  denoted  by  the  velocity  components  U, 

V,  W,  in  the  body  axis  directions;  while  U  ,  V  ,  and  Wa  will  denote  correspond¬ 
ing  velocity  components  with  respect  to  the  (moving)  air  mass.  These  are  related 
to  the  gust  velocity  components  u  ,  v  ,  and  w  seen  by  a  body-axis  observer  as 

O  o  o 

follows: 


U  "  Ug  =  ua 


V  -  v„  =  VQ 
g  a 


W  -  w_  =  W. 
g  a 


(10) 

(ID 

(12) 


The  total  velocity  with  respect  to  the  air  mass  is 


?a  =  (ua2  +  V*  +  "a*) 


1/2 


(13) 


The  true  aerodynamic  angle  of  attack,  a,  and  angle  of  side-slip,  p,  are  defined 
by 


a  =  tan' 


-li 

U 


(14) 


8  =  sin-1  z~ 
V 

a 


(15) 


Other  useful  relations  are 


U_  “  V  cos  8  cos  a 

a  a  K 

V  =  V  sin  8 

a  a  M 

Wa  =  Va  cos  B  sin  a 


(16) 

(17) 

(18) 


In  the  absence  of  atmospheric  turbulence,  the  a  and  p  measurements  corre¬ 
spond  with  the  inertial  values. 


_  +  -1  W 

Gi^  ~~  tan  -q 


(19) 


% 


=  sin 


-II 


v 


4 


(20) 


XB>  yb-  zb  -  THE  body-axis  system  consists  of  right-handed, 

ORTHOGONAL  AXES  WHOSE  ORIGIN  IS  FIXED  AT  THE  NCM> 
INAL  AIRCRAFT  CENTER  OF  GRAVITY.  THE  SYSTEM'S 
ORIENTATION  REMAINS  FIXED  WITH  RESPECT  TO  THE  AIR¬ 
CRAFT,  THE  XB  AND  ZB  AXES  BEING  IN  THE  PLANE  OF 
SYMMETRY.  BECAUSE  THE  EXACT  ALIGNMENT  OF  Xg 
AXIS  IS  ARBITRARY,  HERE  IT  IS  TAKEN  ALONG  THE  BODY 
CENTERLINE  REFERENCE. 


Figure  1 

Aircraft  Axis  System 
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2.  INTEGRATION  OF  THE  AIRSPEED  EQUATIONS 

The  most  direct  method  of  obtaining  is  the  solution  of  three  simultane¬ 
ous,  first-order,  nonlinear  differential  equations. 


U  =  Ax„_  -  QW  +  RV  -  g  sin  6 

'■'to 


V  =  Ay  -  RU  +  PW  +  g  cos  0  sin  0 


W  =  Az._  -  PV  +  QU  +  g  cos  6  cos  0 

eg 


(21) 

(22) 

(23) 


This  requires  measurements  of  Ax  ,  Ay„,  Az„_,  P,  Q,  R,  0,  and  0.  The  equations 

eg  eg  eg 

are  integrated  to  solve  for  U,  V,  and  W,  after  which  Equation  (19)  is  solved. 

One  serious  drawback  to  this  method  is  the  requirement  for  inertial  quality 
computation. 

3.  ALTITUDE  RATE  EQUATION 

Compensated  barometric  altitude  rate,  h^,  assuming  constant  altitude  isobars, 
is  given  by 


b^  =  U  sin  0  -  V  cos  0  sin  0  -  W  cos  0  cos  0 


(24) 


Substituting  inertial  values  for  U,  V,  and  W  similar  to  Equations  (16)  through 
(18)  results  in 


=  V1  cos  0^  cos  sin  0  -  V.^  sin  cos  0  sin  0 

-  \  cos  0,^  sin  cos  0  cos  0 

In  practice,  V&  would  be  substituted  for  V^,  and  would  be  given  by 


(25) 


Vp  p  2  2 

_ cos  0  cos  0  +  sin  6  -  A 

ocj  —  sxn  ~  ppp 

1  cos  0  cos  0  +  sin  0 


where 


ni 

A  =  — - i - +  tan  ^  cos  0  sin  0 

Va  cos  p± 

Equation  (26)  requires  measurement  of  0,  0,  h,  Va,  and 


(26) 


(27) 
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4.  ELIMINATION  OF  U,  U.  V,  V,  W,  AND  W 

Differentiate  Equations  (16)  through  (13)  and  substitute  in  Equations  (1) 
through  (3). 

Va  cos  fa  cos  Va  sin  fa  cos  -  a±  Va  cos  fa  sin 

*  Axcg  -  Q  Va  cos  fa  sin  a,  +R?a  sin  fa  -  g  sin  6  (28) 

Va  sin  fa  +  fa  Va  cos  fa 

*  Aycg  -HVa  cos  fa  cos  a±  +  P  V&  cos  fa  sin  ai 

+  g  cos  0  sin  0  (29) 

Va  cos  fa  sin  ai  -  Va  sin  fa  sin  +  fa  V&  cos  ^  cos  fa 

=  Azcg  -  P  Va  sin  j3i  +  Q  Va  cos  fa  cos  a±  +  g  cos  0  cos  0 

(30) 

Eliminating  P,  Q,  and  R  yields 

—  22  —  2 
Va  cos  fa  cos  -  fa  Va  sin  cos  ^  cos  a± 

•  —  2  —  2  •  — 

-  Qfi  Va  cos  sin  cos  fa  +  V&  sin  fa  +  fa  Va  cos  fa  sin  fa 

—  2  2  •  —  2 

+  Va  cos  fa  sin  a±  -  Va  sin  ^  cos  fa  sin 

♦  —  2 

+  oc±  Va  cos  fa  sin  cos  fa 

«  A*cg  cos  cos  fa  +  Aycg  sin  fa 

+  Azcg  sin  cos  fa  -  g  sin  0  cos  cos  fa 

+  g  cos  0  sin  0  sin  fa  +  g  cos  0  cos  0  sin  cos  j3^  (31) 


4' wl'.  k 


t 


When  collecting  terms,  and  p ^ 

Axcg  cos  a±  cos  /3±  +  AyQg  sin  p±  +  AzQg  sin  a±  cos  p1 

-  /  2  2  .2 
=  Va  (cos  p^  cos  +  sxn  pJ± 

+  cos2  j3^  sin2  +  g^sin  8  cos  cos  p^ 

+  cos  0  sin  0  sin  +  cos  8  cos  0  sin  cos  {}^  (32) 

Equation  (32)  reduces  to  Va,  and  the 


hi 

=  va  +  g  =r 
v. 


(33) 


The  first  term  on  the  right-hand  side  of 
second  term  to  g  Thus 

\  i 

Axcg  cos  ei  cos  pL  +  Aycg  sin  p±  +  AzQg  sin  a±  cos  p±  -  V& 

0  _ 

which  can  be  solved  for  using  measurements  of  Axcg,  Aycg,  Azcg,  p±i  Va,  V 
and  h^. 

5.  EQUATING  LEFT-HAND  SIDES  OF  LIFT  EQUATION  IN  STILL  AIR 
Combining  Equations  (28)  and  (30)  and  reducing, 

(Axcg  +  Va  sin  R  -  g  sin  8)  sin  a.^  -  (Azcg 

-  Va  sin  p^  P  +  g  cos  0  cos  0)  cos  +  V&  cos  p^ 

=  V_  cos  p.  Q  (3*0 

d  -L 

This  can  be  solved  for  a±  in  terms  of  p±,  Axcg,  Azcg,  P,  Q,  R,  V&,  6,  and  0. 


6.  LONGITUDINAL  AND  NORMAL  ACCELERATION  EQUATIONS 

Equations  (28),  the  longitudinal  acceleration  equation,  and  (30),  the 
normal  acceleration  equation,  both  contain  higher-order  terms  and  do  not  appear 
promising  on  an  individual  basis  unless  many  geometric  assumptions  are  made. 

No  steady-state  solution  for  «1  exists,  nor  can  one  be  obtained  by  summing 
these  two  equations. 


7.  PITCHING  MOMENT  EQUATION 

The  aircraft  pitching  moment  equation  is  given  by 


M6  6T  +  a  +  Mq  Q  +  bQ 


=  Q  + 


Y  PQ  +  (p2  “ 

V  y 


(35) 
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This  differential  equation  can  be  solved  for  a  without  explicitly  measuring  Q  if 
M(a)  =  M  a.  The  solution  then  requires  that  P,  Q,  R,  6_,  m,  and  q  be 

0 1  e  i 

measured.  The  method  for  approximate  solution  is  presented  in  Section  III. 

8.  NORMAL  FORCE  EQUATION 

The  apparent  simplicity  of  the  aircraft  normal  force  equation  recommends  its 
use  for  computing  a- 

Azcg  =  Z(“>  +  \  6e  *  \  6T  (36) 


Z(a)  is  a  stability  derivative  that  is  known  with  high  precision.  Its  large 

magnitude  tends  to  relax  the  requirements  for  accuracy  in  representing  Z.  and 

°e 

Z,  .  Measurements  of  q  and  mass  are  also  required. 

°t 

Directly  related  to  the  normal  force  equation  is  the  lift  equation. 


m  (La  “  +  L&e  6e>  =  ^cg  sin  « 


-  Azcg  cos  a  -  Z£ 


cos  a 


(37) 


The  comments  for  the  normal  force  equation  apply  here  also.  Equations  (36)  and 
(37)  are  mainly  different  in  tha^  (37)  requires  an  extra  measurement  (Axcg),  and 
the  solut.-*  on  for  a  Is  made  somewhat  more  complex  because  the  right-hand- side 
terms  are  multiplied  by  functions  of  a. 


9.  SUM  OF  NORMAL  FORCE  AND  PITCHING  MOMENT  EQUATIONS 
Equations  (35)  and  (36)  can  be  summed  to  yield 

Az  -  X  Q  =  Z(a)  +  Z  b  +  M  6™  +  Z  6m  -  X  M*  3  -  X  M(a) 
eg  a'  &ee  6j  i  bj  i  aa  a 

Jx  -  Iz 

-  X„  M  Q  -  X„  M.  +  Xo  -  PQ 

a  q  '  a^e®  a  v 


+  X  (P2  -  R2) 

a  y 


(38) 


If  X  is  chosen  such  that  X  =  (Z.  /M,  ),  then  the  6  term  in  Equation  (38) 


becomes  zero.  Furthermore, 


V  &e 


Ecg  "  K”  «  =  Az 
°e 


(39) 
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where  A’  is  the  output  of  a  normal  accelerometer  located  at  X  =  Z  /M 


V 


and 


Ya  =  Z&  =  0.  Equation  (38)  requires  measurements  of  P,  Q,  R,  A^,  &T,  6T,  and 

mass  along  with  accurate  knowledge  of  Z(a),  Zt,  Z,  ,  M(a),  M,  ,  M*.  ,  Mn, 

o^,7  Og7  7  6g7  a7  q 

Ix’  Iz’  V  8X111  Ixz* 


SECTION  III 


if 


l 


ERROR  ANALYSIS 


Section  II  developed  a  multitude  of  equations  from  which  aircraft  angle  of 
attack  could  be  extracted.  The  following  section  examines  these  equations  with  l 

respect  to  complexity,  practicality,  and  accuracy  in  an  attempt  to  determine  a  1 

feasible  system.  | 

1.  EQUATIONS  YIELDING  INERTIAL  a  j 

Integration  of  Equations  (21)  through  (23)  will  produce  a  high-quality  | 

inertial  a ,  provided  inertial  quality  computer  and  sensors  are  used.  From  a  i 

cost  standpoint  alone,  this  method  is  quite  impractical.  In  addition,  any  1 

attempt  to  simplify  the  equations  will  result  in  cumulative  errors  tv-*>t  quickly  S 

become  intolerable.  I 

B 

Equation  (25),  using  inertial  altitude  rate  and  Euler  angle  measurements, 

appears  to  be  a  promising  algorithm.  It  has  a  steady-state  solution,  thxxs  | 

preventing  cumulative  error  buildup.  Furthermore,  the  sensor  requirements  are  4 

minimal  except  for  a  possible  3  measurement.  This  equation  (repeated  here  for  5 

convenience  with  V„  substituted  for  V.)  is  known  henceforth  as  Method  I.  | 

**  1  1 

hi  =  Va  (cos  f}^  cos  sin  0  -  sin  cos  0  sin  0  | 

J 

-  cos  f}^  sin  a.^  cos  0  cos  0  % 


A  thorough  study  of  this  equation  will  be  presented  following  the  selection  of 
all  other  candidate  methods. 

Equation  (33)  also  appears  to  be  a  useful  algorithm.  It  requires  an  air 
data  system,  three  accelerometers,  and  a  f}  measurement.  Like  Equation  (25),  it 
provides  a  steady-state  solution,  thus  relaxing  the  mechanization  requirements. 
Equation  (33)  will  be  designated  as  Method  II. 

— 

Axcg  cos  cos  +  Aycg  sin  /j±  +  Azcg  sin  oc±  cos  %  =  V&  +  g  — 

Equation  (3*0  has  no  steady-state  solution  for  O)  and.  is  therefore  subject 
to  cumulative  errors.  In  addition,  it  is  considerably  more  complex  than  Methods 
I  and  II,  requiring  accelerometers,  rate  gyros,  Euler  angles,  p,  and  air  data 
measurements.  Thus,  Equation  (3*0  is  not  recommended  for  a  computation. 
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a 


& 


2.  EQUATIONS  YIELDING  TRUE  a 

Equation  (35) }  the  pitching  moment  equation,  is  fairly  complicated  as  given. 
However,  several  assumptions  can  be  made  to  considerably  simplify  the  equation. 
First,  assume  that  P,  Q,  and  R  are  small,  making  their  products  negligible  and 
M(a)  =  a •  Equation  (35)  then  becomes 

bm  +  M  a  +  M.  a  +  Q  +  M.  6  =  Q  (40) 

a  q  &e  e  ' 

Solving  for  a  yields 

« = 4  (k '  \ 61  ■  M«  “  *  \ 15  ■  \  6«)  (4i) 

Taking  some  mathematical  liberties  with  Equation  (4l), 


1/M4  (S  -  Mq)  \/Mi 

S  +  M^  Q  "  S  +  6e 


S  + 


(42) 


indicating  that  Q  need  not  be  measured.  The  equation  does  require  measurements 
of  Q,  &e,  6^,,  mass,  and  Mach  number  and  accurate  knowledge  of  M^,  Mq,  Mq,  M&  , 

and  M&  .  Two  basic  factors  limit  the  usefulness  of  Equation  (42).  First,  and 

most  important,  the  computed  value  of  a  is  an  approximation  at  the  outset. 

Second,  the  equation  requires  accurate  knowledge  of  four  aircraft  stability 
derivatives  as  a  function  of  Mach  and  mass,  which  leads  to  a  complex  mechaniza¬ 
tion.  In  view  of  the  simpler,  more  accurate  methods  available,  Equation  (42)  is 
not  recommended  for  computing  <*• 

The  normal  force  equation,  Equation  (36),  is  given  by 

AZcg  =  ZC«>  +  Z^  6,  +  Z^  6T 

As  was  stated  previously,  the  relative  simplicity  of  this  equation  recommends  its 
use  in  a  computation.  It  is  henceforth  known  as  Method  III  and  will  be  analyzed 
in  greater  depth  later. 

The  lift  equation  is  similar  in  principle  to  the  normal  force  equation  but 
has  additional  complexities.  Since  there  is  nothing  to  recommend  its  use  over 
Method  III,  it  is  discarded.  The  same  can  be  said  for  summing  the  normal  force 
and  pitching  moment  equations. 
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3.  ATMOSPHERIC  TURBULENCE 


Before  examining  the  errors  that  can  develop  due  to  simplifying  assumptions, 
atmospheric  turbulence  effects  should  be  considered.  Both  Methods  I  and  II  com¬ 
pute  inertial  a,  which  is  in  er.’or  due  to  the  contribution  of  gust  velocities  to 
U  and  W_.  Whether  the  c  error  is  serious  depends  on  the  mission  requirements, 
the  severity  of  the  turbulence,  and  the  aircraft  response. 

The  requirements  for  a  high-frequency  a  response  come  primarily  from  the 
weapons  computer  and/or  the  gust  alleviation  system  input.  The  weapons  delivery 
problem  is  complicated  by  the  effect  of  aircraft  velocity.  If  weapons  delivery 
takes  place  at  relatively  high  velocities,  then  effects  of  gusts  on  a  can  be 
sufficiently  small  to  be  ignored.  Whether  inertial  a  is  sufficient  for  weapons 
delivery  will  be  a  function  of  the  individual  aircraft  and  mission  requirements. 

The  requirements  for  a  at  low  speed  are  generally  dictated  by  some  form  of 
speed  command  system  in  which  high-frequency  excitation  is  undesirable .  There¬ 
fore,  a  low  pass  filtered  inertial  a  may  be  adequate  as  an  input.  Once  again, 
the  requirements  of  the  specific  system  must  be  known  to  determine  the  applica¬ 
bility  of  inertial  a  computation. 

To  illustrate  the  problem,  several  simulation  runs  were  cade  using  the  test 


aircraft. 

The 

results  of 

two  such 

i-uns  are 

presented 

below. 

Flight 

Condition 

Mach 

Attitude. 

Turbulence* 

K_ 

"eg 

(BKS 

*11** 

(RMS 

-  degeeg} 

Ei2 

(peak 

ggaass) 

E21  S22 

(RMS  (peak 

degree)  degrees) 

i 

0.25 

Sea 

level 

Slight 

0.1 

0.35 

0.8 

0.34  0.8 

2 

0.9 

Sea 

level 

Very  heavy 

0.35 

0.24 

0.56 

0. 24  0.60 

•Reference  2 

**E. .  and  E, ~  —  Method  I  a  error 

a*: 

Eg-  and  E^  ~  Method  II  a  error 


4.  APPROXIMATIONS 

Of  the  many  different  methods  of  inferential  a  computation,  only  three  remain 
that  are  accurate  and  yet  relatively  simple.  A  disadvantage  (turbulence)  of 
Methods  I  and  II  has  already  been  discussed.  These  two  methods  do.  however,  have 
one  distinct  advantage  over  the  third:  their  universal  applicability  to  any  air¬ 
craft  type  without  calibration.  Method  III,  on  the  ether  hand,  will  require 
accurate  knowledge  of  at  least  two  aerodynamic  coefficients  for  each  aircraft 
type. 
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T..e  approximations  that  can  be  made  are  limited  by  the  wide  range  of  flight 
pp.ai  ■  ^ s.  over  which  a  must  be  computed.  Obviously,  no  small-angle  approxima- 
t.L  'O'  oan  be  made  for  either  0  or  0,  which  can  be  as  large  as  75  degrees. 

The  effects  of  approximations  are  best  illustrated  exercising  the  simu¬ 
lated  aircraft  through  large-angle  maneuvers  wherein  the  principal  equation 
parameters  become  large.  The  aircraft  response  is  shown  in  Figures  2,  3S  and  4. 
As  each  approximation  is  considered,  the  error  produced  should  be  compared  with 
the  aircraft  response  to  gain  insight  into  the  approximation  effects  with  re¬ 
spect  to  each  of  the  equation  parameters. 

Side-slip  angle  can  be  approximated  by  sin  ^  =  0^  and  cos  0^  =  1.  This 
causes  Methods  I  and  II  [Equations  (25)  and  (33)]  to  become,  respectively, 

h.  =  V&  (cos  sin  0  -  0^  cos  0  sin  0  -  sin  cos  0  cos  0)  (43) 

and 

—  ^i 

Axcg  cos  a±  +  0±  Aycg  +  Azcg  sin  =  Va  +  g  --  (44) 

va 

The  effect  of  this  approximation  on  a  computation  is  shown  in  Figure  5*  Since 
the  error  which  occurs  is  small,  this  approximation  is  useful  for  all  but  the 
most  stringent  system  requirements. 
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Figure  3 

Aircraft  9  and  0  Response 
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Aircraft  and  k  Response 
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Figure  5 

Errors  Due  to  Side-Slip  Angle  Approximation 


If  side-slip  angle  is  completely  ignored,  Equations  (25)  and  (33)  simplify 
to 

h.  =  Va  (cos  sin  8  -  sin  a ^  cos  8  oos  0)  (45) 

and 

•“  hi 

Axcg  cos  a±  +  Azcg  sin  a±  =  V&  +  g  —  (46) 

Y& 

This  last  assumption  results  in  considerable  error  under  certain  conditions. 
The  error  occurs  in  Method  I  when  is  large  (10  degrees)  and  a  roll  angle  is 
present.  As  much  as  1.5  degrees  of  error  can  result  if  0  is  only  5  degrees. 
Larger  bank  angles  (30  to  40  degrees)  can  produce  considerably  large  errors  (8  to 
10  degrees).  This  is  presented  graphically  in  Figure  6. 

Method  II  will  be  in  error  if  0  is  large  and  a  lateral  acceleration  is 
present.  For  the  particular  case  shown  in  Figure  6,  the  error  in  Method  II  is 
smaller  than  that  of  Method  I,  but  this  is  not  always  true.  Larger  lateral 
accelerations  will  produce  equivalent  errors  in  Method  II. 

Geometric  approximations  for  sin  a  and  cos  a  simplify  the  mechanization  re¬ 
quirements  considerably.  However,  the  errors  produced  are  unacceptable  for  all 
but  the  lowest  quality  system.  These  errors  are  illustrated  in  Figures  7>  8,  9, 
and  10. 


Method  III  does  not  lend  itself  to  geometric  simplifications.  Those  param¬ 
eters  which  are  not  critical  include  C_  and  Z  thrust.  Depending  on  aircraft 

6e 

type,  Z  thrust  may  be  zero.  The  representation  of  C  in  the  mechanization  may 


require  only  a  simple  constant  gain.  On  the  other  hand,  Cz(a)  must  be  precisely 
represented,  including  the  effects  of  any  high-lift  devices. 
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Figure  6 

Errors  Due  to  Ignoring  Side-Slip 
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Error  in  Method  I 
when  Sin  a  =  a  and  Cos  a  =  1 
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Figure  10 
Error  in  Method  II 
when  Sin  a  ~  a  and  Cos  ct  =  1 
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SECTION  r/ 


SENSOR  ANALYSIS 


1.  INTRODUCTION 

To  gain  insight  into  the  effects  of  sensor-induced  errors,  the  study  has 
been  separated  into  maneuvering  (dynamic)  and  non-maneuvering  (static)  flight. 
The  static  case  uses  reduced  equations  to  reveal  the  basic  accuracy  to  be  ex¬ 
pected  from  each  computing  method.  The  dynamic  equations,  which  illustrate  the 
effects  of  each  equation  parameter,  nave  been  evaluated  on  a  worst-case  basis 
since  a  completely  analytical  approach  proved  ■•menlightenlng. 

Ideally,  the  sensor  inputs  to  the  a  computer  should  have  zero  error.  Any 
improvement  in  basic  sensor  accuracy  will  Improve  the  a  computation.  However, 
it  is  beyond  the  scope  of  this  program  to  Improve  basic  sensor  designs  so  as  to 
increase  their  accuracy.  We  can  at  best  point  out  those  sensors  which  are 
critical  to  the  accurate  computation  of  a.  Operations  involving  sensor  outputs 
such  as  filtering,  quickening,  and  mode  conversion  are  considered  where  improve¬ 
ments  In  a  computation  are  obtained. 

2.  STATIC  ERRORS 

Examining  the  first  two  methods  of  a  computation  under  conditions  of  wings- 
level,  constant-altitude  flight  reveals  their  basis  for  computing  steady-state 
a.  In  the  absence  of  Input  offsets  in  static  flight,  Equations  (25)  and  (36) 
reduce  to,  respectively, 


-sin  cos  6  +  cos  a^  sin  6=0  (47) 

and 

Ax„„  cos  a.  +  Az  sin  a.  =  0  (48) 

eg  i  eg  I 

Equation  (47)  merely  states  that  for  vings-level,  constant-altitude  flight, 
a.;  =  6.  Thus,  the  static  computation  of  a^  in  Method  I  Is  simply  the  output  of 
the  vertical  reference.  This  implies  that  the  c  computation  can  be  no  better 
than  the  vertical  reference  accuracy  when  using  Method  I. 

Equation  (48)  has  reduced  in  Method  II  to  the  arctangent  of  the  ratio  of 

Ax  to  Az  .  But  in  static  flight,  Ax  =  -g  sin  8  and  Az  =  g  cos  8,  and 
eg  eg  eg  eg 

again,  a.^  =  8.  In  this  case,  the  accuracy  of  the  a  computation  Is  dependent  cn 
the  basic  accelerometer  accuracies  and  their  alignments  with  the  X  and  Z  body 
axes. 
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The  third  computation  method  for  a  does  not  reduce  in  the  static  case,  but 
the  equation  allows  a  mathematical  approach  to  the  error  analysis,  which  will  be 
discussed  subsequently. 

Table  I  shows  the  effects  of  static  sensor  errors  on  the  three  methods  of  a 
computation  for  the  six  flight  conditions.  In  each  case,  one  sensor  wa3  given  a 
5-percent  offset  from  the  trim  value,  while  the  other  sensor  outputs  remained 
accurate.  The  values  presented  reveal  that  in  Methods  I  and  II  the  a  computa¬ 
tion  accuracy  is  directly  proportional  to  the  accuracy  of  the  input  sensors.  In 
Method  III  it  is  apparent  that  three  of  the  seven  parameters  are  not  critical; 
those  being  thrust,  C  ,  and  elevator  position.  This  is  of  particular  impor- 

6e 

tance  with  respect  to  thrust  since  it  may  be  difficult  to  accurately  measure  in 
complex  engines. 
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TABLE  I 

ANGLE-OF-ATTACK  ERROR  AT  TRIM  DUE  TO  STATIC  SENSOR  ERRORS 


» 


A 


3.  DYNAMIC  ERRORS 

The  analysis  of  dynamic  errors  considers  the  unreduced  Equations  (25),  (33), 
and  (36)  in  the  presence  of  maneuvering  flight.  A  worst-case  study  of  all  three 
methods  is  followed  by  a  mathematical  error  analysis  of  Method  III. 

3.1  Worst-Case  Study 

A  worst-case  analysis  of  sensor  errors  has  been  used  due  to  the  relatively 
fruitless  results  obtained  from  a  purely  mathematical  approach.  Some  special 
situations  in  which  the  equations  reduce  to  a  reasonable  complexity  will  be 
mathematically  analyzed  later. 


As  an  example  of  the  complexity  resulting  from  a  mathematical  approach,  con¬ 
sider  the  solution  of  Equation  (25)  using  the  identity 


cos 


/  2  a/2 

=  I  1  -  sin  J 


Equation  (25)  then  becomes 


+  sin  p^  cos  0  sin  0 


-sin  cos  0  cos  0  +  ( 1  -  sin 


^1  -  sin2  a±)‘ 


sin  0 


M) 

cos 

(50) 


When  solved  for  0^,  this  yields 


a.  =  sin 


■I  A  cos  0  cos  0  +  sin  0  (cos2  0  cos  2  0  +  sin2  0 


cos2  0  cos2  0  +  sin2  0 


l2) 


1/2 


(5D 


where 


A  = 


tan  p^  cos  0  sin  0 


(52) 


The  sensitivity  of  to  the  individual  parameters  is  determined  by  taking  the 
partial  derivative  of  Equation  (5l)  with  respect  to  the  subject  parameter.  The 
derivative  of  arc sin  is  given  by 


/  2»"l/2 

da  =  (1  -  u  /  du 


(53) 


Equation  (53)  applied  to  Equation  (51)  for  any  of  the  parameters  leads  to  a 
highly  complex  expression  that  provides  no  insight  into  the  parameter  effects. 
Furthermore,  if  a  quantitative  evaluation  is  attempted,  realistic  values  must  be 
assigned  to  each  parameter.  This  either  limits  the  analysis  to  one  or  two 
points,  or  inflates  the  task  beyond  reasonable  bounds. 
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The  worst-case  analysis  procedure  consisted  of  examining  the  effects  of  an 
erroneous  sensor  output  during  large  aircraft  maneuvers.  Flight  Conditions  (FC) 
1  and  3  were  used  to  provide  both  low-  and  high-performance  cases.  Large  bank 
angles  and  side-slip  angles  were  produced  using  FC-1,  whereas  FC-3  provided 
large  normal  and  lateral  accelerations.  In  each  instance,  the  maneuvers  were 
limited  to  produce  realistic  parameter  values.  "Worst  case"  must  be  qualified 
by  stating  that  although  the  maneuvers  may  not  represent  a  perfect  worst-case 
situation  with  respect  to  a  computation,  they  do  exercise  all  aircraft  param¬ 
eters  beyond  normal  conditions  without  resorting  to  an  unrecoverable  attitude. 

In  each  simulation  run,  one  sensor  was  given  a  5-percent  scale  error  while 
all  others  remained  accurate.  The  influence  of  sensor  inaccuracy  was  then  ob¬ 
served  by  recording  the  error  between  the  computed  and  true  a.  The  aircraft 
response  for  FC-1  (Mach  0.25  at  sea  level)  is  shown  in  Figure  11.  The  errors 
induced  in  each  of  the  three  computing  methods  are  presented  in  Figures  12 
through  14.  The  high-performance  case,  FC-3  (Mach  0.9  at  sea  level),  is  shown 
in  Figures  15  through  18. 
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Figure  12 
Sensor  Errors 

Flight  Condition  1,  Method  I 
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Figure  13 
Sensor  Errors 

Flight  Condition  1,  Method  II 


36 


(DEGREES) 


Figure  13  -  Continued 


TIME  (SECONDS) 


TIME  (SECONDS) 


TIME  (SECONDS) 


TIME  (SECONDS) 


Figure  15 

Aircraft  Response,  Flight  Condition  3 
(a)  a  and  B 
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Figure  18 
Sensor  Errors 

Flight  Condition  3»  Method  III 
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Figure  18  -  Concluded 


p  Is  the  only  parameter  in  Method  I  that  does  not  require  high  accuracy. 

Bank  angle  is  not  critical  at  high  speeds  but  does  require  accuracy  at  low 
speeds.  If  any  parameter  can  be  considered  the  most  critical,  it  is  the  ver¬ 
tical  reference  in  the  longitudinal  axis,  6.  For  the  major  portion  of  most 
flights,  a  ~  0;  and  if  8  is  not  accurate,  a  will  seldom  be  correct.  Altitude 
rate  and  airspeed  are  important  primarily  when  there  is  a  negative  flight  path 

angle  and  6  is  close  to  aero.  In  this  case,  y~~ZT  and  a  ~  v*  view  of  the 

Va 

magnitude  of  error  caused  by  each  of  the  sensors,  at  least  1-percent  accuracy  Is 
required  for  all  but  the  side-slip  angle  measurement.  A  p  measurement  of  5  per¬ 
cent  should  be  sufficient  in  any  case. 

Two  parameters  in  Method  II,  Ay  and  p,  require  only  5-percent  accuracy. 

eg 

The  remaining  sensors,  Axcg,  AzCg’  hi»  Va’  and  va}  re<lu:I-re  a  1-percent 
measurement. 

•  _ 

The  large  errors  occurring  in  Method  II  with  errors  In  h^,  Vg,  \T&,  Aycg,  and 

Ax  are  correlated  with  normal  acceleration  (Az  )  passing  through  or  close  to 
c^.r  eg 

;erc.  The  discontinuities  are  explained  mathematically  by  examining  Equation 
(S')  when  Azcg  is  0.  Sol.'Ing  for  a^, 


(54) 


where 


g^i  i 

K  =  —  +  v„  -  sin  p.  Ay 


i  "'eg 


The  sensitivity  to  errors  in  K  is  found  by  taking  the  partial  derivative  of 


When  Az„  is  close  to  zero,  a*  is  small,  and  from  Equation  (54),  Ax^  -  must 

Cg  J- 

be  small.  This  makes  the  denominator  of  Equation  (56)  small,  and  the  sensitiv¬ 
ity  of  to  changes  in  K  extremely  large.  The  solution  to  this  problem  Is  to 
inhibit  computation  of  a.,  when  using  Method  II  If  Az  is  close  to  zero.  This 
is  discussed  in  Section  V. 


3.2  Method  III  Sensitivity  Analysis 

The  independence  of  the  parameters  in  Method  III  permits  a  useful  mathemati¬ 
cal  sensitivity  analysis.  The  particular  item  of  interest  is  the  RMS  error  that 
would  occur  if  the  sensors  were  each  in  error  by  a  specified  amount.  The  RMS 
error  is  given  by 


Method  III  has  been  evaluated  at  trim  for  FC-1  and  FC-3  assuming  each  of  the 
sensors  and  parameters  is  in  error  by  1  percent.  The  resulting  computed  angle- 
of-attack  errors  were  0.20  degree  for  a  =  11.2  degrees,  and  0.034  degree  for 
a  -  2.1  degrees.  A  reasonable  limit  to  the  accuracy  of  Method  III  seems  to  be 
approximately  2  percent,  assuming  that  all  parameters  can  be  measured  to  1 
percent. 

4.  SENSOR  IMPROVEMENTS 

The  following  discussion  compares  sensor  requirements  with  respect  to  a  com¬ 
putation  against  state-of-the-art  sensor  performance  characteristics.  It  is 
already  apparent  from  the  static  error  and  sensitivity  analyses  that  the  0.1- 
degree  accuracy  requirement  will  be  difficult  to  meet  with  any  of  the  three  com¬ 
putation  methods.  For  this  reason,  a  new  set  of  constraints  Is  proposed  based 
on  the  a  measurement  requirements  for  the  latest  generation  of  high-performance 
aircraft.  These  requirements  are  listed  below. 

Accuracy  Threshold  Time  Constant 
Aircraft  Velocity  (degree)  (degree)  ( second) 

90  to  125  knots  ±0.5  0.2  0.075 

125  knots  to  Mach  3-0  ±0.2  0.1  0.075 

4.1  Vertical  Gyro 

The  vertical  gyro  provides  measurement  of  Euler  angles  9  and  0.  Except  in  a 
low-grade,  low-cost  system,  it  is  inadequate  for  angle-of-attack  computation  due 
to  its  precession  under  sustained  accelerations. 

The  long-term  precession  can  be  overcome  by  combining  Methods  I  ar.d  II. 
Method  I  will  provide  good,  short-term  accuracy  (using  a  vertical  gyro)  while 
Method  II  is  best  in  the  long  term.  Thus,  a  mechanization  of  both  methods  — 


* 


passing  the  first  through  a  high  pass  filter  and  the  second  through  a  low  pass 
filter  and  then  summing  the  two  —  will  provide  a  high-quality  inertial  a  with 
good  dynamic  response.  Low  passing  Method  II  has  the  added  advantage  of  filter¬ 
ing  the  V  signal,  which  tends  to  he  noisy.  Figure  19  is  a  block  diagram  of  a 
simplified  mechanization  of  this  combination. 


4.2  Inertial  Navigation  System 

Euler  angles  8  and  0  can  be  obtained  from  an  Inertial  Navigation  System 
(INS)  with  sufficient  accuracy  for  the  a  computation.  It  is  expected  that  the 
availability  of  an  INS  will  increase  with  each  new  generation  of  aircraft.  The 
INS  does  not  suffer  from  long-term  accelerations  and  is  therefore  preferred  over 
the  vertical  gyro  In  Method  I. 


4.3  Air  Data  System 

The  air  data  system  is  required  for  all  three  a  computation  methods.  The 
necessary  outputs  are  Ya^r,  Vair,  Mach,  and  p.  Th<-»  air  data  system  accuracy 
and  frequency  response  are  marginal  with  respect  to  a  computation.  Particularly 
at  the  lower  speeds,  the  airspeed  inacc  racy  can  be  significant. 


At  altitude,  the  lag  in  the  true  airspeed  measurement  can  be  considerable  If 
the  pressure  transducers  are  remote  from  the  source.  It  would  be  desirable  to 
compensate  for  this  lag.  Let  the  measurement  of  true  airspeed  bs  Vffi  where 


Vm  = 


T(p)  s  +  1 


(58) 


The  time  constant,  T,  will  be  principally  a  function  of  density  for  a  given  in¬ 
stallation.  The  measurement  can  be  partially  compensated  for  plumbing  lag, 
introducing  a  V  term  as  shown  below. 

a 


V 


a 


v  T(P)  s  +  1  +  Tl  Kg  -  g  5ln  6) 

m  T-j_  s  -r  1  T-^  s  +  1 


(59) 


Airspeed,  altitude  rate,  Mach  number,  and  air  density  are  standard  air  data 
outputs.  Airspeed  rate  is  obtained  by  pseudo  differentiation  of  the  airspeed 
signa] . 

The  accuracy  of  Methods  I  and  II  could  be  greatly  improved  in  turbulence  if 
air-mass-related  altitude  rate  were  available.  Unfortunately,  for  all  other 
applications  inertial  altitude  rate  is  the  required  quantity,  and  the  Air  Data 
System  (ADS)  is  designed  to  minimize  the  effects  of  turbulence  on  h^.  Enhancing 
these  turbulence  effects  would  require  relocating  the  static  pressure  part  to  a 
turbulence- sensitive  position. 
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Dynamic  pressure,  ■*  p  V  ,  can  be  calculated  within  the  a  computer  if  it 
is  not  available  as  an  ADS  output. 

4.4  Side-Slip  Angle 

The  measurement  of  p  presents  an  interesting  problem  in  that  a  probe  or  vane 
must  be  located  in  the  airstream.  If  [5  is  not  available,  Methods  I  and  II  can 
exhibit  more  than  1.5  degrees  of  error  during  conditions  of  large  p  and  small 
bank  angles  or  lateral  accelerations. 

Side-slip  can  be  calculated  internally  in  a  manner  similar  to  a  computation. 
The  least  complex  method  is  the  Y-force  Equation  (4). 


tycg  -  ^ 


3  +  Y. 


The  relative  simplicity  of  this  equation  would  recommend  its  use.  However,  It 
has  one  serious  flaw:  the  value  of  Ya  may  be  extremely  difficult  to  obtain  and 
is  generally  known  with  poor  precision.  Nevertheless,  the  equation  still  re¬ 
mains  the  best  choice  if  p  is  to  be  calculated. 


4.5  Accelerometers 

Method  II  requires  accelerometers  located  at  the  aircraft  center  of  gravity 
in  all  three  body  axes.  Most  types  of  accelerometers  will  ra^et  the  requirements 
of  a  computation.  One  item  that  deserves  particular  attention  is  the  misalign¬ 
ment  between  the  accelerometer’s  sensitive  axis  and  the  desired  body  axis.  The 
error  can  be  geometrically  compensated  if  the  angles  of  misalignment  are  known. 
(This  subject  is  well  documented  In  Reference  3.)  Additionally,  the  accelerom¬ 
eters  must  be  mounted  as  close  to  the  aircraft  center  of  gravity  as  is  practical. 


4.6  Elevator  Position 

Measurement  of  aircraft  elevator  position  does  not  present  any  formidable 
problems.  The  measurement  Is  generally  available  with  more  than  sufficient 
accuracy  on  most  advanced  aircraft. 


4.7  Thrust 

Measurement  of  thrust  is  required  only  if  the  thrust  line  is  not  parallel  to 
the  X  body  axis.  Even  slight  misalignments  may  be  tolerable  since  the  measure¬ 
ment  is  not  critical.  In  the  event  that  thrust  misalignment  is  significant, 
throttle  position  may  provide  sufficient  accuracy  for  a  computation. 
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4.8  Mass 

No  direct  measurement  of  mass  exists  other  than  takeoff  mass  on  at  least  one 
of  the  latest  generation  of  cargo  aircraft.  Mass  can  be  calculated  using  fuel 
flow  measurements  or  the  lift  equation.  Use  of  the  lift  equation  appears  to  be 
the  most  promising  technique  and  will  be  discussed  in  Section  V.  In  any  system 
using  calculated  mass,  a  takeoff  initial  condition  may  have  to  be  inputed 
manually  or  through  a  strain  gauge/landing  gear  measurement  system. 

5.  QUICKENING  DEVICES 

The  term  "quickening  devices"  has  been  used  in  conjunction  with  angle-of- 
attack  systems  development  by  the  U.S.  Navy.  The  basic  purpose  of  the  quicken¬ 
ing  device  was  to  minimize  the  lag  in  the  pilot  response  during  carrier  landings. 
This  was  accomplished  by  putting  a  lead  network  in  the  display  to  reflect  the 
result  of  the  pilot  input  before  the  aircraft  fully  responded.  Quoting  from  the 
flight  evaluation  final  report*:  "The  purpose  of  the  concept  was  to  enable  the 
pilot  to  control  a  (or  analogously  airspeed)  during  une  carrier  approach  by 
minimizing  the  inherent  time  lag  in  the  pilot-airplane-.:  control  loop.  The 
modification  was  termed  ’instantaneous  angle-of-attack  display  (IAAD) 1  with  a 
’quickening’  characteristic  of  the  display." 

As  it  turned  out,  quickening  devices  did  not  appear  to  be  a  promising  Inno¬ 
vation.  Again  quoting  from  the  flight  test  report:  "The  IAAD  is  unsuitable  for 
service  use  from  the  standpoint  of  monitoring  ease  and  resultant  pilot  perfor¬ 
mance  because  of  three  basic  discrepancies,  i.e.,  excessive  system  sensitivity, 
lack  of  situation  display  and  inadequate  system  Input  design  for  airplane  re¬ 
sponse  characteristics."  Further  difficulties  of  quickened  a  were  detailed  in 
the  flight  test  report. 

Quickening  devices  as  they  relate  to  angle-of-attack  computation  do  not 
appear  to  offer  any  advantages  with  respect  to  system  performance. 

6.  SENSOR  STUDY  CONCLUSIONS 

In  general,  airframe  angle  of  attack  can  be  accurately  computed  using  cur¬ 
rently  available  on-board  sensors.  Since  the  weakest  link  in  the  computing 
system  is  the  ADS,  any  efforts  to  improve  sensor  operation  with  respect  to  a 
computation  should  be  concentrated  on  air  data  functions. 


*Reference  4. 
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Filtering  except  for  producing  airspeed  rate  is  confined  to  low  pass  net¬ 
works  designed  to  lioit  sensor  frequency  response  and  noise  output.  The  air¬ 
speed  rate  signal  is  derived  froa  the  airspeed  signal  through  a  lead-lag  network. 

Ultimately,  the  sensor  requirement  will  be  dictated  by  the  mission  require¬ 
ments  with  respect  to  a.  Systems  not  requiring  a  for  weapons  computation  or 
gust  alleviation  may  considerably  relax  the  need  for  high-quality  sensors. 
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SECTION  V 

MECHANIZATION  STUDIES 


1.  DISCUSSION  OF  THE  EQUATIONS 

1.1  Method  I,  the  Euler  Angle  Equation 

The  Euler  angle  equation,  (25),  has  a  serious  difficulty  which  limits  its 
usefulness  in  high-quality,  angle-of-attack  computation.  For  accurate  computa¬ 
tion,  a  knowledge  of  side-slip  angle  (p)  is  required.  Since  one  of  the  goals  of 
an  angle-of-attack  computer  is  to  avoid  the  need  for  an  a  vane,  it  is  unrealis¬ 
tic  to  assume  that  a  vane  is  available  for  p  measurement.  Therefore,  p  must  be 
computed,  and  an  accurate  computation  of  p  is  at  least  as  difficult  as  computing 
a.  A  possibility  might  be  to  use  Method  II,  the  accelerometer  equation,  to  com¬ 
pute  side-slip  angle  and  Method  I  to  compute  angle  of  attack  (or  vice  versa). 

On  examination,  however,  both  Methods  I  and  II  reduce,  during  straight  and  level 
flight,  to  the  trivial  case 

0  sin  p^  +  0  cos  p^  =  0  (60) 

The  most  practical  method  available  for  computation  of  p  is  the  Y-force 
equation,  (4-).  This  equation  requires  a  knowledge  of  mass,  the  computation  of 
which  is  quite  complex. 

It  is  concluded  that  Method  I  is  most  useful  in  a  simplified  system  where 
side-slip  angle  can  be  neglected.  If  it  is  assumed  that  is  small  enough  to 
permit  small-angle  approximations,  the  following  equation  results. 


cos  6  cos  0 


(sine't) 


Since  the  limits  on  8  and  0  are  ±75  degrees,  there  is  no  problem  of  division  by 
zero.  The  airspeed,  V  ,  is  zero  only  when  the  airplane  is  on  the  ground  and  the 
system  is  not  operational.  A  block  diagram  (Figure  20)  shows  one  possible 
mechanization. 

1.2  Method  II,  the  Accelerometer  Equation 

Some  of  the  comments  for  Method  I  also  apply  to  Method  II,  Equation  (33). 

The  computed  angle  of  attack  is  inertial  rather  than  air-mass  related,  and  the 
need  exists  for  computation  of  p-  In  addition,  airspeed  rate,  which  is  not 
normally  available  from  the  air  data  computer,  is  required.  Although  it  is  . 
noisy  and  somewhat  difficult  to  use,  pseudo  rate  can  be  derived  and  used  for  V&. 
True  airspeed,  Va,  would  be  used  in  mechanization  rather  than  inertial  velocity. 
V,  because  of  availability. 


For  these  reasons,  Method  II  is  not  recommended  for  high-quality  a  computa¬ 
tion.  However,  in  a  simpler,  more  economical  system,  a  simplified  version  of 
this  method  can  be  used.  By  assuming  Aycg  sin  =  0,  cos  p^  =  1,  sin  =  a^, 
and  cos  =  1,  the  following  equation  is  obtained. 

= + g  t  -  <62> 

This  equation  can  be  solved  explicitly  for  aj_. 

A  difficulty  occurs  when  attempting  to  divide  by  zero.  Az.  can  go  to  zero 

eg 

in  flight  during  a  zero  g  maneuver.  This  can  be  handled  by  logically  detecting 

when  a  low  value  of  Az„„  exists  and  clamping  the  angle-of -attack  output  to  the 

eg 

value  previously  computed,  or  possibly,  by  clamping  to  the  known  zero-lift  angle 
of  attack  for  that  flight  condition. 

Static  errors  (errors  during  trimmed  flight)  are  produced  by  the  approxi¬ 
mations  cos  =  1  and  sin  If  this  approximation  is  not  made,  the  equa¬ 

tion  cannot  be  solved  explicitly  for  but  must  be  solved  implicitly. 

Implicit  computation  is  a  closed-loop  technique  whereby  the  computed  variable 
converges  to  the  proper  value  through  a  continuous  adjustment.  A  common  example 
is  the  performance  of  division  by  using  a  multiplier  in  a  high-gain  feedback 
loop.  As  with  any  closed-loop  system,  stability  can  be  a  problem.  In  the 
divider  example,  if  the  divisor  changes  sign,  static  instability  results.  An 
approach  to  implicit  computation  designed  to  minimize  some  error  criterion  (for 
example,  mean  square  error)  can  avoid  the  stability  problem  and  minimize  errors 
at  difficult  points  such  as  the  indeterminate  zero  divided  by  zero. 

Assume  that  the  given  equation  takes  the  form 

A(t)  -  B(os)  =  0  (63) 

where  it  is  desired  to  evaluate  a«  In  the  implicit  solution,  a  will  be  adjusted 
to  minimize  the  mean  square  error  between  A(t)  and  B(a).  In  other  words, 

A(t)  -  B(a)  =  e  (64) 

the  error  to  be  minimized.  Define  the  mean  square  error 

E  =  e2  (65) 
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Minimum  E  occurs  when  the  derivative  of  E  with  respect  to  a  equals  zero.  Con¬ 
sidering  Figure  21,  which  shows  E  and  3E/3a  plotted  as  a  function  of  a  in  the 
vicinity  of  minimum  E,  note  that 

•  When  3E/3a  is  negative,  a  is  too  small 

•  When  3E/6a  is  positive,  a  is  too  large 

•  When  3E/3a  =0,  a  is  optimum 

It  follows  that  a  will  seek  its  optimum  value  provided 


da  v-  ,3E 
dt  "*■  3a 


where  K  is 
Taking 


some  positive  constant. 

the  derivative  Equation  (65)} 


3E  _ 

3a 


(66) 


(67) 


and  from  Equation  (64) 


3e  _  _  8B(oO 
3a  3a 


Thus, 


i  =  acef^ 


(68) 

(69) 


An  example  of  this  technique  is  shown  in  Figure  22,  an  implementation  of 
Method  II  assuming  ^  =  0.  Writing  the  equation  in  the  form  of  Equation  (63), 


& 


V  +  V 
a 


♦.-( 


Azcg  sin  a±  +  AxQg  cos  a± 


)-° 


For  this, 


S|Sl=Aj 

3a 


■cg  00  3  ai  -  ^cg  sin  ai 


(70) 


(7D 


The  equation  mechanized  in  Figure  22  is 


m.  =  K  (Azcg  cos  -  Axcg  sin 


■■)W 


#  +  V  - 


Va  -  Azcg  sin  «i 


-  Ax„„  cos 
eg 


0 


“i)  (72) 


This  technique  provides  a  stable  computation  when  Azcg  passes  through  zero 

and  reverses  polarity.  Although  the  accuracy  of  solution  is  sensitive  to  sensor 

errors  at  Az  =  0,  the  computational  method  is  valid, 
eg 
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Figure  21 

Variation  of  E  and  ^  as  a  Function  of 

Oft 
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Figure  22 

Implicit  Computation,  Method  II 


1.3  Method  III,  the  Z-Force  Equation 

Equation  (36),  the  z-force  equation,  is  the  best  method  for  high-quality, 
angle-of -attack  computation.  As  for  mechanization,  difficulty  arises  in  that 
aircraft  mass  must  be  continuously  computed,  and  aerodynamic  coefficients,  which 
vary  with  Mach  and  airplane  configuration,  must  be  computed. 

The  most  practical  method  for  computing  mass  is  through  the  z-force  equation. 
However,  this  equation  cannot  be  used  to  compute  both  mass  and  angle  of  attack. 
The  lift  equation  could  be  used  to  compute  mass,  but  this  is  just  the  z-force 
equation  rotated  by  the  angle  a,  and  is  not,  for  practical  purposes,  sufficiently 
independent  of  the  z-force  equation  to  permit  derivation  of  two  unknown 
quantities. 

The  solution  proposed  is  to  use  Method  I  or  II  to  compute  inertial  angle  of 
attack,  and  use  this  in  the  z-force  equation  to  compute  mass.  Then  with  the 
z-force  equation,  a  good,  high-quality,  dynamic  angle  of  attack  can  be  computed. 
Statically,  the  resulting  angle  of  attack  can  be  no  better  than  the  previously 
computed,  independent  angle  of  attack  used  in  computing  mass.  Therefore,  this 
must  be  quite  accurate  statically. 

With  Method  II  used  to  compute  <x1}  the  simplifying  assumptions 
Aycg  sin  f}^  =  0  and  cos  0^=1  are  permissible  because  they  are  dynamic  terms. 
However,  sin  and  cos  must  be  computed,  or  at  least  more  closely  approxi¬ 
mated  than  simply  sin  and  cos  =  1.  As  shown  in  subsequent  paragraphs, 

if  a  Digital  Differential  Analyzer  (DDA)  is  used  for  the  computations,  deriva¬ 
tion  of  the  sine  and  cosine  is  a  simple  matter.  But  If  a  special  purpose, 

whole -number  computer  or  an  analog  computer  is  used,  the  approximations 
2  3 

cos  a  =  1  -  a  /2  and  sin  a  -  a  -  a/6  are  easier  to  mechanize. 

Figure  23  shows  the  block  diagram  of  this  approach.  The  complexity  is  such 
that  it  is  feasible  to  consider  special  purpose  digital  computing  methods  to 
mechanize  the  system.  Indeed,  the  accuracy  required  to  obtain  a  maximum  error 
of  0.1  degree  when  a  ranges  from  0  to  20  degrees  is  better  than  0.5  percent. 

This  is  not  achievable  using  analog  circuitry  at  reasonable  cost. 

Computation  of  the  aerodynamic  coefficient,  Cz(a,  Ma),  can  be  achieved  by 
tables  and  interpolation  in  a  digital  computer.  Alternatively,  if  a  functional 
relationship  can  be  empirically  determined,  expressing  C„  as  a  function  of 
powers  of.  a  and  Mach,  the  function  can  be  computed  directly.  Asstuning  this 
function  can  be  expressed  with  sufficient  accuracy  and  without  undue  complexity, 
this  approach  would  be  well  suited  for  either  DDA  or  whole-number  digital 
computation. 
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As  previously  mentioned,  care  must  be  taken  when  dividing  by  quantities  that 
can  go  to  zero.  In  this  system,  only  Azcg  can  go  to  zero  when  the  system  is 
operating.  This  can  best  be  handled  logically  by  detecting  small  or  zero  Azcg 
and  holding  the  appropriate  variables  constant. 

The  thrust  term  must  be  included  if  the  thrust  line  is  not  the  same  as  the 
x  body  axis.  It  is  assumed  that  thrust  is  a  function  of  throttle  position  only. 
A  simple  polynomial  of  second  or  third  degree  should  be  suffi Jiently  accurate 
since  the  thrust  term  is  small,  and  the  resulting  error  in  computed  a  will  be 
insignificant. 


MACH,  FLAPS,  WING  SWEEP,  ETC. 
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Figure  23 
High- Accuracy, 
Angle- of- Attack  Computer 
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2.  ANALOG  COMPUTATION 

Analog  techniques  provide  the  best  means  for  producing  an  economical  angle- 
of-attack  system.  This  is  due  to  the  high  cost  of  the  analog-digital  interface 
equipment  as  veil  as  the  greater  complexity  of  digital  systems.  Analog  cir¬ 
cuitry  has  accuracy  limitations  that  make  it  unsuitable  for  the  high-quality 
system,  but  it  is  entirely  adequate  for  applications  such  as  transports  and 
cargo  airplanes. 

The  simplified  versions  of  Methods  I  or  II  could  readily  be  mechanized 
using  analog  circuitry,  with  certain  advantages  over  digital  hardware.  The 
almost  infinite  resolution  of  the  analog  voltages  permits  the  use  of  extremely 
high  gains  in  implicit  computation  loops.  In  digital  implicit  loops,  finite 
resolution  of  the  data  results  in  limit  cycle  oscillations  at  high-loop  gains. 
Another  advantage  of  analog  systems  is  the  ease  of  interface:  analog  sensors 
and  displays  are  widely  used  at  present,  while  their  digital  counterparts  are 
just  beginning  to  enter  service.  Analog  systems  will  continue  to  hold  an 
economic  advantage  at  the  interface  until  the  majority  of  the  interfacing  units 
are  digital. 

The  multiplier  is  the  key  component  in  the  analog  circuitry.  The  multi¬ 
pliers  found  most  in  practice  are  pulse- width  modulation  and  quarter- square 
types,  although  the  variable  transconductance  type  is  becoming  more  widely  used. 
The  quarter- square  multipliers  are  standard  in  analog  computer  facilities  be¬ 
cause  of  their  high-accuracy  capabilities,  but  they  are  very  expensive.  Pulse- 
width  multipliers  are  used  extensively  in  flight  control  equipment  for  gain 
changing.  They  are  relatively  inexpensive,  but  accuracy  is  limited  to  approxi¬ 
mately  1  percent.  They  also  have,  the  disadvantage  of  low  bandwidth  since  the 
output  must  be  filtered  to  remove  the  carrier  frequency  component. 

Both  hybrid  and  microcircuit  multipliers  of  the  variable  transconductance 
type  have  been  recently  introduced  on  the  market.  These  multipliers  only 
achieve  accuracies  of  about  1  percent  but  they  are  very  economical,  and  it  is 
reasonable  to  assume  that  their  accuracy  will  improve  as  the  market  for  multi¬ 
pliers  expands.  Several  major  semiconductor  companies  are  developing  low-cost 
multipliers  of  this  type. 

The  accuracy  of  multipliers  over  a  temperature  range  is  a  much  more  dif¬ 
ficult  problem  than  that  of  operational  amplifiers.  In  an  amplifier,  stable 
resistors  are  used  to  define  the  gain;  but  in  a  multiplier,  gain  (scale  factor) 
is  at  least  partially  determined  by  the  stability  of  semiconductor  components. 
In  addition,  there  are  offset  errors  that  vary  with  temperature.  Due  to  the 


lack  of  detailed  specifications  on  recently  available  models,  it  is  difficult 
to  determine  actual  performance  over  temperature  without  testing. 

A  circuit  diagram  of  Method  II  simplified  (/S^  =  0,  sin  a ^  =  a cos  a-j_  =  1) 
is  shown  in  Figure  24.  Equation  (62)  is  multiplied  through  by  V_  and  Az„_  and 
solved  implicitly,  resulting  in  a  slightly  simpler  mechanization  than  solving 
for  a^  explicitly.  Explicit  computation  requires  two  divide  operations,  which 
are  really  implicit  computations  using  multipliers.  Hence,  the  technique 
employed  is  simpler  because  it  has  only  one  implicit  loop  rather  than  two. 

The  problem  of  static  instability  due  to  polarity  changes  of  Az  is 

eg 

handled  by  switching  polarity  in  the  loop.  Two  polarity  detectors  are  used 
having  a  region  of  overlap  around  Az  =  0.  Thus,  in  this  region  the  input  to 

C  0 

the  integrator  is  shorted  to  ground,  and  the  output  remains  constant  during  the 
zero  g  maneuver. 

This  mechanization  requires  eight  encapsulated  circuit  modules,  including 
two  multiplier  modules  and  a  power  supply  module.  The  type  of  multiplier  is 
not  specified,  but  accuracy  is  specified  to  be  1  percent  over  the  temperature 
range,  without  external  calibration.  It  is  assumed  that  such  a  multiplier  will 
be  available  in  the  near  future  at  a  cost  of  $100  each  in  quantity. 

The  packaged  angle-of-attack  computer  using  the  analog  mechanization  of 
Method  II  would  weigh  approximately  1.0  pound,  occupy  a  volume  of  32  cubic 
inches,  and  consume  3  watts  of  power.  The  sales  price  is  estimated  to  be  $600 
to  $800  each. 

Method  I  is  slightly  more  complex  than  Method  II.  It  requires  Scott-Tee 
transformers  to  convert  synchro  pitch  and  roll  data  to  sines  and  cosines,  and 
good  demodulators  to  convert  the  sines  and  cosines  to  dc  voltages.  Three 
multipliers  are  required  instead  of  two.  Twelve  modules  are  required;  they 
weigh  approximately  1.5  pounds,  occupy  a  volume  of  48  cubic  Inches,  and  consume 
4  watts.  The  cost  is  estimated  to  be  $800  to  $1000  each. 

3.  DIGITAL  DIFFERENTIAL  ANALYZER 

When  mechanizing  a  problem  that  is  small  or  moderate  in  complexity  and 
requires  accurate  nonlinear  or  functional  computations,  the  DDA  is  a  prime  con¬ 
tender.  A  DDA  is  a  special  purpose  digital  computer  that  performs  incremental 
computation  by  simply  updating  the  values  previously  computed  rather  than  re¬ 
computing  the  entire  function  each  cycle.  Any  operation  that  can  be  expressed 
in  terms  of  differential  equations  can  be  readily  mechanized. 
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In  many  respects,  the  DDA  is  like  an  analog  computer.  The  basic  computing 
element  is  an  integrator.  A  DDA  system  can  be  described  by  means  of  a  "patch" 
diagram  showing  the  interconnection  of  computing  elements  to  solve  a  given  set 
of  differential  equations.  The  information  flow  between  integrators  is  in  the 
form  of  pulses  representing  changes  to  be  made  to  the  problem  variables. 

3.1  Basic  DDA  Operation 

The  symbol  for  the  integrator  is  shown  in  Figure  25.  The  basic  integrator 
consists  of  two  registers  and  two  adder/subtracters.  The  Y-register  contains 
the  current  value  of  the  dependent  variable  (Y),  and  the  R-register  holds  the 
remainder  resulting  from  the  previous  operation.  The  inputs  and  outputs  of  the 
integrators  are  incremental  variables  but  in  the  figure  are  written  as  the 
differential  variables  which  they  approximate.  When  a  dy  increment  appears  at 
the  input,  Y  is  incremented  by  a  constant  scale  factor  q  (decremented  if  the 
sign  of  cy  is  negative).  When  a  dx  increment  appears,  Y  is  added  to  R  (sub¬ 
tracted  if  the  sign  of  dx  is  negative).  Whenever  the  value  of  R  +  Y  exceeds  p 
(the  capacity  of  the  R-register),  an  output  dz  pulse  is  generated  and  the 
R-register  holds  the  remainder.  When  the  size  of  the  increments  is  small  com¬ 
pared  to  the  whole-word  problem  variables,  the  numerical  solution  of  differen¬ 
tial  equations  can  be  quite  accurate. 

A  DDA  servo  is  used  to  sum  incremental  signals  and  provide  the  output  in 
incremental  form.  It  is  also  used  as  a  high-gain  element  in  the  same  manner 
that  an  open-loop  amplifier  is  used  in  an  analog  computer.  In  operation,  the 
contents  of  a  register  in  the  DDA  servo  is  incremented  by  each  of  several  dy 
inputs.  A  zero  detection  circuit  examines  the  contents  of  the  register  and 
provides  an  output  whenever  the  contents  is  anything  other  than  zero.  Sometimes 
the  output  is  gated  and  the  polarity  controlled  by  a  dx  input  signal.  The  servo 
symbol  is  similar  to  the  integrator  symbol  with  the  word  "Servo"  written  in  it. 

3.2  Reversibility 

Due  to  the  finite  levels  of  quantization,  a  DDA  computing  system  will 
exhibit  systematic  error  when  the  polarity  of  the  independent  variable  is 
reversed.  The  error  may  build  up  indefinitely  during  computation  and  result 
in  large  accumulative  errors.  Fortunately,  there  are  methods  of  handling  this 
problem,  making  it  possible  to  eliminate  the  problem  completely  in  most  cases 
and  reducing  it  to  very  small  levels  in  others.  In  general,  a  system  of  DDA 
integrators  will  be  reversible  without  error  if  each  DDA  unit  is  reversible 
and  if  all  computations  in  the  DDA  system  due  to  a  single  input  pulse  are 
completed  before  a  second  input  pulse  arrives.  These  conditions  are  always 
sufficient  but  may  not  be  necessary  for  a  particular  system. 
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Y  =  Y0  +  qy  =  Yq  +  q  J*  dy 

WHERE 

p  =  OVERFLOW  POINT 
q  =  Y  INCREMENT 
Y0  =  INITIAL  VALUE  OF  Y 


Figure  25 

Symbol  and  Equation  for  DDA  Integrator 
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The  best  way  to  make  a  DDA  unit  reversible  is  to  mechanize  it  to  perform 
trapezoidal  integration  rather  than  rectangular  integration.  This  has  the  added 
advantage  that  more  accurate  integration  is  performed.  However,  true  trape¬ 
zoidal  integration  is  not  always  possible  because  it  can  only  occur  when  there 
is  no  time  delay  associated  with  computation  of  the  inputs  to  the  integrator. 

In  any  closed- loop  computation,  there  will  be  at  least  one  integrator  whose  in¬ 
puts  were  computed  during  the  previous  computation  cycle  and  hence,  have  a  time 
delay  of  one  sampling  period.  Even  so,  a  form  of  trapezoidal  integration  that 
greatly  reduces  any  errors  can  be  mechanized. 

Controlled  timing  satisfies  the  requirements  that  all  computations  due  to 
a  single  input  pulse  be  completed  before  a  second  input  pulse  arrives.  If  the 
computer  is  fast  enough  to  complete  all  computations  within  the  time  between  two 
successive  increments  of  any  variable  changing  at  its  maximum  rate,  the  require¬ 
ment  is  met.  This  is  sometimes  very  difficult  to  accomplish  in  practice. 

Achieving  reversible  computation  depends  to  a  great  extent  on  the  type  of 
function  being  computed.  Of  the  functions  most  used  in  the  alpha  computer,  the 
sin/cos  function,  polynomial  evaluation,  and  multiplication  functions  are  easily 
made  reversible,  while  division  is  very  difficult. 

3.3  Incremental  Computation 

There  are  certain  consequences  of  performing  incremental  computation  that 
must  be  considered.  First,  in  order  to  get  the  computation  started,  all  the 
data  registers  must  be  filled  with  their  proper  initial  conditions.  This  means 
that  fixed  constant  numbers  are  stored  in  the  machine  so  that  when  power  is 
initially  applied  these  numbers  are  automatically  set  into  the  registers. 

Second,  the  sensor  input  data  to  the  machine  must  be  converted,  not  only 
from  Analog- to-Digital  (A/D)  but  also  to  incremental  form.  This  is  commonly 
done  in  conjunction  with  a  standard  A/D  converter.  The  increment  generator 
(Figure  26)  consists  of  storage  space  for  each  input  variable  in  a  data  memory, 
an  adder/subtracter,  and  a  comparator.  An  analog  input  to  the  A/D  converter  is 
compared  to  the  accumulated  value  of  the  increments  of  that  quantity  held  in 
storage.  If  they  are  not  equal  an  increment  is  generated,  positive  If  the  in¬ 
put  is  larger  and  negative  if  the  input  is  smaller  than  the  stored  value.  The 
increment  is  then  added  to  the  stored  value,  which  is  put  back  in  storage  to 
await  the  next  conversion  cycle. 

If  an  error  is  made  in  a  DDA  due  to  noise,  power  supply  transient,  etc, 
the  DDA  will  be  in  error  for  all  subsequent  time  due  to  its  incremental  nature. 
Subsequent  computations  are  merely  changes  to  what  has  gone  on  before,  and  when 
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an  error  goes  in,  it  remains  in.  There  are  two  ways  to  attack  this  problem. 

One  is  to  detect  the  existence  of  certain  check  points  of  an  input  and  set  the 
computer  directly  to  the  corresponding  values.  For  example,  if  sine  and  cosine 
of  an  input  angle  are  being  computed,  the  zero  degree  point  might  be  designated 
the  check  point.  Whenever  the  angle  passes  through  zero,  the  sine  can  be  set  to 
zero  and  the  cosine  to  one  immediately.  The  incremental  computation  then 
carries  on  from  that  point.  If  any  errors  have  been  made,  they  are  corrected; 
if  not,  no  change  is  made  in  the  data. 

Another  method  of  correcting  any  transient  errors  that  may  have  occurred 
is  to  set  the  computer  registers,  including  the  increment  generator  storage 
registers,  to  a  specific  set  of  points  with  the  outputs  and  all  intermediate 
variables  correctly  corresponding  to  the  values  set  in  the  increment  generator 
registers.  Then  the  computer  slews  up  to  the  set  of  values  corresponding  to 
the  actual  inputs  at  the  maximum  computer  slew  rate.  This  operation  can  be 
repeated  periodically  to  ensure  that  no  error  remains  in  the  computer  longer 
than  a  fixed  amount  of  time. 

These  two  techniques  require  additional  hardware  over  that  needed  to  per¬ 
form  the  actual  computations.  However,  they  have  both  been  successfully  used 
in  DDA  practice. 

3. V  Hate  Limiting 

The  frequency  response  of  a  DDA  is  limited  in  that  the  problem  variable 
can  change  by  only  one  increment  per  computing  cycle.  If  the  input  variables 
are  changing  at  a  very  high  rate  or  at  such  a  rate  that  the  computed  output 
must  vary  rapidly,  the  DDA  may  contribute  a  considerable  lag.  The  two  factors 
that  affect  the  maximum  rate  the  computer  can  handle  are  the  computer  cycle 
time  and  the  size  of  the  increments. 

Since  a  problem  variable  may  change  by  one  increment  per  cycle,  the  larger 
the  Increment,  the  greater  the  rate  of  change  of  the  variable  for  a  given  cycle 
time.  However,  large  increments  cause  poor  computational  accuracy  and  lead  to 
limit  cycle  oscillations  in  closed-loop  systems. 

The  computer  cycle  time  depends  on  the  computer  organization  and  on  the 
clock  rate  at  which  the  circuitry  operates.  Assuming  the  system  is  operating 
at  the  maximum  clock  frequency  that  the  circuitry  can  tolerate,  the  frequency 
response  will  depend  only  on  the  organization. 

Organization  refers  to  the  degree  of  parallelism  in  the  computer. 
Arithmetic  processing  can  be  either  serial  or  parallel,  and  processing  of  the 
integrators  either  sequential  or  simultaneous.  Economy  demands  a  minimum  of 


80 


i  libm,  but  this  is  achieved  at  the  expense  of  speed»  Most  DDA's  use 

;’erial  arithmetic,  but  processing  of  the  integrators  varies  from  all-sequential 

all-parallel.  If  a  configuration  of  32  integrators  is  not  fast  enough  for  a 
particular  application,  the  DDA  can  be  organized  into  two  groups  of  16  inte¬ 
grators  operating  in  parallel  to  double  the  speed.  The  cost  would  actually  be 
quite  small  since  the  memory  requirements  are  the  same  for  both  systems.  Only 
an  additional  arithmetic  unit,  consisting  of  several  adders  and  logical  circuits, 
would  br  required. 

3.5  DDA  Meehan ’ zation 

The  DDA  mechanization  of  the  system  in  Figure  23  is  developed  for  esti¬ 
mating  cost,  size,  and  performance.  The  design  is  based  on  the  use  of  Metal- 

On-Silicon  (MOS)  shift  registers  for  the  data  registers,  and  read-only  memory 
for  program  control  and  constant  storage.  Hence,  the  arithmetic  is  serial.  It 
is  assumed  that  the  integrators  will  be  processed  sequentially. 

A  clock  rate  of  5  MHz  will  be  used.  A  word  length  of  20  bits  is  dictated 
by  resolution  requirements  (10  data  bits,  1  sign  bit,  and  9  bits  for  scale 
factor  resolution).  The  resulting  time  required  to  perform  one  integration  is 
4  microseconds. 

Figure  27  is  a  DDA  patch  diagram  of  the  computer.  As  indicated  by  the 
figure,  30  integrators  and  13  DDA  servos  are  needed.  A  machine  built  with  only 
the  minimum  number  of  units  would  be  totally  inflexible  and  without  growth 
capability*  Hence,  a  minimum  of  32  integrators  and  16  servos  is  used  for  the 
estimates  in  this  study. 

Given  the  clock  rate,  number  of  integrators,  and  word  length,  the  slewing 
rate  capability  of  the  machine  can  be  determined.  For  an  all- sequential  machine, 
the  computer  cycle  time  is  4  microseconds  per  integrator  multiplied  by  48  inte¬ 
grators,  or  192  microseconds.  Assuming  an  increment  of  0.025  degree,  the  out¬ 
put  slewing  rate  is  then  130  degrees  per  second.  Thus,  it  appears  that  a 
sequential  process  DDA  will  handle  the  ar.gle-of-attack  problem. 

3.6  System  Organization 

Figure  28  shows  the  organization  of  the  DDA.  It  consists  mainly  of  inter¬ 
face  equipment,  the  arithmetic  unit,  timing  and  control  circuits,  and  various 
memories. 

The  interface  equipment  contains  a  multiplexer  capable  of  handling  a  mini¬ 
mum  of  eight  inputs,  an  A/D  converter,  an  increment  generator,  and  a  D/A 
converter.  The  A/D  converter  can  be  time  shared  for  D  '  .  conversions  if  used 
with  sample- and- hold  output  amplifiers. 
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Figure  27 
DDA  Patch  Diagram 


83 


SCALE  FACTOR 


SERVO 


VARIABLE 

SCALE 

FACTORS 


F -719-2-41  (2) 


Figure  27  -  Concluded 
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Figure  28 
DDA  Block  Diagram 
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The  arithmetic  unit  consists  of  four  serial  adder/subtracters ,  overflow 
detection  logic,  and  zero  detection  logic.  This  circuitry  is  time  shared  by  all 
integrators  and  servos. 

The  timing  circuitry  consists  of  two  counters  (a  bit  counter  and  a  word 
counter)  and  two  decoders.  The  bit  counter  and  decoder  define  bit  times  within 
a  word,  while  the  word  counter  and  decoder  define  word  times  within  a  computer 
cycle. 

There  are  five  distinct  memories  in  the  system.  The  MOS  data  memory  is  a 
pair  of  long  shift  registers  containing  all  the  Y  and  R  registers  of  the  inte¬ 
grators  and  servos.  The  increment  generator,  whole-word  memory  is  a  long  shift 
register  storing  the  whole-word  digital  inputs  from  the  A/D  converter.  The 
increment  memory  stores  the  outputs  from  the  increment  generator  and  each  DDA 
integrator  and  servo.  All  of  the  memories  discussed  so  far  are  shift  register 
types.  The  program  memory  and  scale  constant  memory  are  solid-state,  read-only 
types.  The  program  memory  provides  addresses  to  the  increment  memory  at  the 
proper  times  to  select  outputs  from  integrators  as  inputs  to  other  integrators 
required  to  execute  the  program.  Scale  constants  are  also  gated  by  timing 
signals  into  the  arithmetic  unit  as  needed. 

The  physical  characteristics  of  the  DDA  are  estimated  as  follows: 

e  Size  175  cubic  inches 

•  Weight  6  pounds 

•  Power  25  watts 

The  cost  estimate  for  the  DDA  is  based  on  the  use  of  large-scale  MOS  mem¬ 
ories,  and  Transistor- Transistor  Logic  (TTL)  including  Medium-Scale  Integrated 
(MSI)  circuits.  The  parts  are  full  military  temperature  range  types,  and  the 
prices  assume  quantity  purchase.  On  this  basis,  the  sales  price  of  the  DDA 
system  in  production  would  be  approximately  $4000. 
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4.  WHOLE-NUMBER  COMPUTER 

A  small,  special  purpose,  whole-number  digital  computer  would  be  ideal  for 
performing  the  angle-of-attack  computations.  The  only  question  concerns  the 
economics  of  this  approach  when  compared  with  the  DDA,  which  is  normally  consid¬ 
ered  the  simplest  digital  approach  to  a  small  problem.  A  detailed  analysis  of 
the  computer  requirements  has  been  made,  permitting  realistic  estimates  of  the 
hardware  requirements.  The  following  paragraphs  describe  a  hypothetical  com¬ 
puter  that  has  the  capability  required  for  computing  angle  of  attack  and  some 
reserve  capacity  for  self-test  routines  and  growth.  This  model  is  used  for  es¬ 
timating  purposes  comparing  it  to  the  DDA  design  of  Paragraph  3» 

4.1  Advantages  of  Whole-Number  Computation 

Whole-number  computers  have  several  important  advantages  over  incremental 
computers.  A  whole-number  computer  has  greater  flexibility  to  perform  a  wide 
variety  of  computations,  logical  operations,  and  self-test  routines.  The  basic 
DDA  solves  only  differential  equations,  and  must  have  considerable  extra  hard¬ 
ware  to  provide  any  different  type  of  operation. 

The  analog-digital  interface  is  simpler  for  the  whole-number  computer 
because  no  increment  generator  is  required.  In  addition,  since  each  output  is 
recomputed  each  cycle,  any  error  caused  by  noise  or  power  supply  transient  will 
be  corrected  the  next  cycle. 

The  whole-number  computer  will  normally  require  much  less  data  memory  than 
the  DDA.  In  the  latter  unit,  all  problem  variables  are  stored  in  data  memory  so 
they  can  be  updated  each  computer  cycle.  In  a  whole-number  computer,  most  com¬ 
puted  variables  are  required  only  temporarily;  that  is,  until  they  are  used  in 
some  other  computation.  Then  they  are  discarded,  and  the  memory  locations  in 
which  they  were  stored  are  available  for  new  variables. 

There  is  no  reversibility  problem  with  a  whole-number  computer,  nor  is  there 
rate  limiting,  except  where  a  rate  limit  is  deliberately  programmed.  The  fre¬ 
quency  response  is  limited  only  by  the  computation  cycle  rate  and  sampling 
theory  considerations. 

4.2  Computational  Requirements 

Referring  to  Figure  23,  the  majority  of  computations  are  multiplications  and 
additions.  There  are  also  a  sizable  number  of  scaling  operations  required;  that 
is,  multiplication  by  constants.  Scaling  can  often  be  accomplished  by  shift 
operations,  which  save  a  great  deal  of  time.  A  few  divisions  are  required,  and 
if  the  hardware  is  included  to  provide  for  multiplication,  the  mechanization  of 
the  divide  instruction  costs  very  little  more. 
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The  z -force  coefficient  can  be  computed  by  polynomial  approximation  as  in 
the  DDA,  or  alternatively,  by  using  storage  tables  and  interpolating  to  deter¬ 
mine  the  value  between  stored  points.  Interpolation  is  a  common  technique  for 
function  evaluation,  with  several  approaches  that  can  be  used.  The  interpola¬ 
tion  routine  consists  of  passing  a  polynomial  through  several  points  surround¬ 
ing  the  input  point.  The  polynomial  may  be  of  any  degree.  The  first-degree 
(linear)  interpolation  is  easiest  to  compute  but  requires  the  largest  number  of 
stored  data  points.  Higher-degree  interpolations  can  achieve  the  required 
accuracy  with  fewer  stored  points,  but  the  computation  is  more  complex  and  re¬ 
quires  more  stored  program  instructions. 

The  limiting  case  is  where  there  are  no  stored  data  points  at  all,  only  co¬ 
efficients  of  a  high-degree  polynomial  describing  the  entire  curve.  This  is 
the  polynomial  approximation  approach  adopted  in  the  DDA  system.  The  polynomial 
approximation  is  better  than  tables/interpolation  in  the  whole-number  computer 
as  well  as  the  DDA  because  it  is  very  easy  to  program  the  polynominal  evaluation 
for  a  reasonably  smooth  curve  such  as  C„(a).  If  the  function  had  a  large  number 

£j 

of  wiggles  and  sharp  bends,  then  a  tables/interpolation  approach  would  be  more 
efficient. 

4.3  Computer  Organization  (Figure  29) 

A  simple,  straightforward  computer  is  envisioned  to  take  maximum  advantage 
of  newly  emerging  solid-state  memory  technology.  The  computer  is  a  serial, 
binary,  2» s  complement  machine  operating  at  a  basic  clock  rate  of  2  MHz.  It  can 
perform  an  addition  in  16  microseconds  and  a  multiplication  in  144  microseconds, 
including  access  to  memory.  There  are  10  Instructions,  some  of  which  can  be 
modified  by  use  of  a  hardware  index  register.  The  characterisitics  of  the  com¬ 
puter  are  summarized  in  Table  II. 

It  is  estimated  that  the  computer  can  perform  the  angle- of- at tack  program  at 
a  rate  xceeding  50  times  per  second,  with  a  considerable  margin  for  growth. 

This  is  based  on  the  use  of  serial  arithmetic,  16-bit  words,  and  a  2-MHz  clock 
rate.  The  total  number  of  instructions  to  be  executed  is  estimated  at  150, 
broken  down  as  follows: 

Additions*  120  at  16  microseconds  =  1,920 

Multiplications  85  at  144  microseconds  =  12,240 

Divisions  5  at  288  microseconds  =  1,440 

*Add,  Subtract,  Load,  and  Store  15,600  microseconds 

15,600  microseconds  per  cycle  is  equivalent  to  about  64  iterations  per  second. 
The  computer  should  not  be  made  faster  than  necessary  since  this  will  add  to  its 
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TABLE  II 

WHOLE-NUMBER  COMPUTER  CHARACTERISTICS 


Data  Type 


Number  of  Instructions 
Computing  Time 

Memory 

Input/Output 
Hardware  Type 

Weight 

Size 

Power 


Serial 

Binary 

2*  s  Complement 
Fixed  Point 
Ten 

Add:  16  Microseconds 

Multiply:  144  Microseconds 

256  Words,  10- Bit,  Read-Only 

31  Wbrds,  16-Bit,  Scratchpad 

Direct  Access  to  Scratchpad 

MOS  Memories 

TTL  Logic 

4.5  pounds 

1L;  Cubic  Inches 

20  Watts 


4.3*1  Instruction  Format 


An  instruction  word  consists  of  10  bits,  subdivided  as  shown  below: 


10 

Index 


8  7 

Order 


4  3 

Address 


Bits  1  through  5  specify  the  location  in  data  memory  to  which  the  instruction 
refers;  bits  6  through  9  identify  the  order  to  be  performed.  A  "one"  in  bit 
position  10  causes  the  contents  of  the  index  register  to  be  added  to  the  address 
before  executing  the  instruction.  If  bit  10  is  "zero",  the  address  is  not 
modified.  A  list  of  instructions  is  presented  in  Table  III. 


TABLE  III 
INSTRUCTION  LIST 

Bit  No.  Instruction 

1  Load 

2  Store 

3  Add 

4  Subtract 

5  Multiply 

6  Divide 

7  Jump 

8  Jump  on  Negative 

9  Shift  Right  (sign  extend,  end  off) 

10  Shift  Left  (circularly) 
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Figure  29 

Whole-Number  Computer 
Block  Diagram 
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*+»3.2  Control  Unit 

The  control  unit  obtains  the  instruction  to  be  executed  from  storage  and 
then  controls  the  various  transfer  gates  necessary  to  execute  this  instruction. 
The  unit  consists  of  an  index  register,  an  index  adder,  a  program  counter,  an 
instruction  register,  timing  counters,  and  a  sequencer.  The  sequencer  is  the 
most  complex  part  of  the  control  unit 5  it  decodes  the  instruction  and  generates 
the  proper  sequence  of  timing  signals  to  execute  the  instruction.  An  index 
register,  which  requires  very  little  hardware,  is  included  to  facilitate  the 
setting  up  of  loops.  This  is  highly  useful  in  programming  polynomials  and  in 
self- test  routines. 

The  data  memory  is  small  and  can  be  addressed  with  only  5  bits.  The  instruc¬ 
tion  word  is  10  bits  in  length:  5  for  operation  (op)  code  and  index,  and  5  for 
addressing.  There  are  256  program  memory  words,  which  takes  a  minimum  of  8  bits 
for  addressing.  One  of  the  32  data  memory  addresses  is  reserved  for  use  when 
addressing  the  program  memory.  If  this  address  appears  with  an  op  code,  the 
computer  takes  the  next  word  in  the  instruction  sequence  as  the  operand  address. 
Thus,  a  10-bit  address  is  available  for  addressing  the  program  memory. 

A  real-time  clock  provides  a  signal  that  interrupts  the  program  and  makes 
it  return  to  the  beginning  of  the  computation  cycle.  This  internally  generated 
signal  occurs  once  every  20  milliseconds.  The  interrupt  operation  not  only 
serves  to  time  the  required  problem  solution  rate,  but  also  provides  a  auto¬ 
matic  recovery  from  improper  program  loops  which  may  arise  due  to  a  temp  rary 
internal  failure. 

4.3.3  Arithmetic  Unit 

The  arithmetic  unit  performs  the  various  operations  required  in  the  solution 
of  a  given  problem.  Two  arithmetic  shift  registers  and  a  serial  adder/ subtracter 
are  the  essential  components.  Addition  is  accomplished  by  adding  a  number  from 
the  data  memory  to  the  number  in  the  accumulator.  This  sum  is  held  in  the 
accumulator. 

When  a  multiplication  is  performed,  the  number  in  memory  is  the  multiplicand, 
and  the  number  in  the  accumulator  is  transferred  to  the  M-Q  register  and  used  as 
the  multiplier.  The  multiplicand  is  successively  added  to  the  accumulator  con¬ 
tents  and  shifted  in  accordance  with  the  "one"  bit  in  the  multiplier.  Upon  com¬ 
pletion  of  the  operation,  the  most  significant  half  of  the  product  is  held  in 
the  accumulator. 

When  a  divide  operation  is  performed,  the  contents  of  the  accumulator  is  the 
dividend  and  the  specified  word  from  memory  is  the  divisor.  Division  is  accom¬ 
plished  using  the  nonrestoring  method.  When  the  operation  is  complete,  the 
quotient  is  in  the  accumulator. 
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4-. 3. 4  Memory 

The  memory  is  divided  into  two  parts:  the  data  memory  and  the  program  mem¬ 
ory.  Both  are  solid-state  types,  well  suited  to  large- scale- integration 
techniques. 

A  small,  random-access  data  memory  is  used  as  a  scratchpad.  It  includes 
certain  addressable  arithmetic  registers,  temporary  storage  for  intermediate 
computations,  and  buffer  storage  for  the  input/output  unit.  It  contains  31 
sixteen-bit  words. 

The  program  memory  is  read-only  storing  the  fixed  program  and  problem  con¬ 
stants.  It  consists  of  256  ten-bit  words. 

4.3.5  Input/Output  Unit 

A  straightforward,  flexible  input/output  unit  has  space  assigned  in  the  data 
memory  location  for  each  A/D  converter  input  word  and  each  D/A  output  word.  The 
Input/Output  (I/O)  unit  accesses  the  data  memory  directly,  putting  in  words  as 
they  are  converted  and  taking  out  data  as  needed  for  the  D/A  converter.  The 
program  uses  the  data  in  these  specified  locations  without  reference  to  the  I/O 
unit.  The  data  transfers  occur  in  synchronism  with  the  computer  clock,  but  the 
converters  themselves  operate  asynchronously.  They  need  operate  only  fast 
enough  to  make  each  conversion  approximately  once  per  computation  cycle. 

4.4  Cost  Estimate 

The  cost  estimate  is  based  on  the  use  of  MOS  solid-state  memories,  TTL  logic 
including  MSI,  full  military  temperature -range  parts,  and  quantity  purchase 
prices.  The  current  sales  price  of  the  digital  system  in  production  would  be 
approximately  $3000. 

The  failure  rate  of  conventional  transducers  deserves  farther  discussion. 
Vanes  and  probes  typically  advertise  MTBF's  .on  the  order  of  10,000  hours.  These 
numbers  never  reflect  what  really  happens  to  a  probe  or  vane  after  mounting  on 
an  aircraft.  The  transducer  is  used  as  a  handle,  step,  brace,  etc.  No  mounting 
position  has  yet  been  devised  to  eliminate  the  handling  and  environmental  prob¬ 
lems.  Thus,  it  would  appear  that  the  10,000-hour  calculated  MTBF  must  be  judged 
with  extreme  prudence. 

Semiconductor  prices  have  shown  a  considerable  shift  downward  each  year.  An 
accurate  prediction  of  price  in  3  to  4  years  is  difficult  to  obtain,  but  it 
would  not  be  unreasonable  to  reduce  the  $3000  figure  by  an  amount  ranging  from 
$250  to  $500.  At  the  same  time,  the  recent  sophisticated  vane  and  probe  prices 
have  approached  $2000.  With  the  more  stringent  environmental  requirements  ex¬ 
pected  for  future  generation  aircraft,  vane  and  probe  prices  may  well  approach 
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the  computer  cost.  When  the  additional  aspect  of  reliability  is  injected,  the 
computer  (MTBF  10,000  hours)  cost  effectiveness  compares  favorably  with  con¬ 
ventional  transducers  (MTBF  ranging  from  37  to  2500  hours). 

One  of  the  largest  items  of  expense  is  the  analog-digital  interface  equip¬ 
ment.  Eventually  this  can  be  greatly  reduced  by  interfacing  with  only  digital 
equipment.  Digital  air  data  computers  and  accelerometers  having  pulse  outputs 
presently  exist.  In  addition,  instruments  and  displays  that  accept  digital  in¬ 
puts  have  been  developed.  If  only  a  digital  interface  is  required,  the  price 
estimate  for  the  angle-of-attack  computer  is  $2200. 

5.  SAMPLED-DATA  COMPUTER 

In  Reference  5j  a  sampled-data  computer  technique  is  described  in  which  the 
primary  problem  variables  are  represented  as  analog  voltages,  and  amplifier 
gains  are  represented  by  digital  and  discrete  control  functions.  The  potential 
advantages  of  this  technique  are  increased  flexibility  and  programmability  of  a 
digital  computer  wi  hout  the  penalty  of  A/D  and  D/A  conversion  equipment.  The 
digital  gain  control  function  is  achieved  by  converting  the  digital  input  into 
a  pulse  width  that  is  used  to  pulse-width  modulate  the  analog  signal.  Multipli¬ 
cations  are  performed  by  averaging  the  pulse-width  modulated  signal  in  an  inte¬ 
grator  circuit.  The  digital  pulse-width  modulator  is  time  shared  by  all  the 
using  analog  amplifiers  in  the  system. 

Nonlinear  functions  can  be  computed  in  the  sampled-data  computer.  By  incor¬ 
porating  an  analog  comparator  in  the  basic  computer,  it  is  possible  to  evaluate 
integrals  of  the  form 


where  x  and  w  are  independent  analog  variables,  and  y(x' )  is  any  function  of  the 
machine  variable  x* .  With  such  integrals,  it  is  possible  to  perform  a  great 
variety  of  functions  such  as  logarithms,  exponentials,  and  sine-cosines,  as  well 
as  analog  multiplication  and  division. 

The  sampled-data  computer  is  not  suitable  for  angle-of-attack  computation 
because 

•  It  is  not  very  accurate  when  compared  to  a  good-quality  analog  system;  it 
is  limited  in  accuracy  to  approximately  5  percent. 

•  The  digital  pulse-width  modulator,  even  when  multiplexed  and  time  shared, 
adds  greatly  to  the  complexity.  The  only  advantage  is  flexibility  of  pro¬ 
gramming,  which  has  very  little  value  in  either  the  simple-  or  high- 
quality  angle-of-attack  computer  systems. 
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SECTION  VI 


CONCLUSIONS 


There  are  three  practical  methods  of  computing  aircraft  angle  of  attack 
from  internally  mounted,  previously  available  sensors.  The  first  method  uses 
a  vertical  gyro,  an  air  data  computer,  and  a  side-slip  angle  transducer  to  pro¬ 
duce  inertial  a.  The  second  method  requires  three  body-axis  mounted  accelerom¬ 
eters,  an  air  data  computer,  and  a  side-slip  angle  transducer  and,  like  the 
first  method,  yields  inertial  a.  The  third  method  uses  a  normal  accelerometer, 
an  elevator  position  transducer,  a  flap  transducer,  a  throttle  position  trans¬ 
ducer,  and  an  air  data  computer  to  produce  true  a.  A  calculation  of  mass  along 
with  accurate  knowledge  of  two  key  airframe  parameters  is  also  required. 

Of  these  three  approaches,  only  Method  III  has  the  ultimate  capability  of 
computing  errorless  angle  of  attack  under  all  conditions  of  flight;  the  first 
two  methods  exhibit  errors  when  flying  in  turbulence.  On  the  other  hand, 

Method  III  is  considerably  more  complex  than  the  other  two,  requiring  knowledge 
of  aircraft  parameters,  a  mass  calculation,  and  a  sophisticated  mechanization. 
The  two  inertial  methods  are  easily  mechanized  and  are  universally  adaptable 
to  any  aircraft  type. 

It  is  concluded  that  a  combination  of  Methods  I  and  II  would  be  adequate 
for  transport  applications.  Fighter-type  aircraft  would  require  a  mechanization 
of  the  third  method,  while  bomber  types  fall  in  between  but  lean  to  the  third 
method. 

The  initial  cost  of  an  angle-of-attack  computer  will  probably  exceed  that 
of  a  vane  or  probe  transducer.  However,  the  superior  reliability  of  the  com¬ 
puter  is  expected  to  more  than  compensate  for  the  price  differential  when 
maintenance  costs  are  considered.  The  angle-of-attack  computer  also  offers 
potentially  higher  accuracy  than  the  external  devices. 
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APPENDIX 


AIRCRAFT  DATA 


The  data  presented  in  this  appendix  represents  the  aerodynamic 
characteristics  of  a  high-performance  tactical  fighter  aircraft.  It  is  a 
composite  from  several  different  data  sources  and  does  not  necessarily  represent 
a  particular  aircraft  model.  Of  the  eleven  flight  conditions  available,  six 
were  used  in  the  study  program  and  only  those  six  are  presented. 
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TABLE  IV 

AERODYNAMIC  DATA  FOR  SIMULATED  AIRCRAFT 


Paraaatar 

Plight  Condition 

1 

3 

7 

8 

10 

u  1 

h(ft) 

0 

0 

15,000 

35,000 

0 

35,000 

Mach 

0.25 

0.9 

1.1 

0.6 

0.154 

2.5 

P( slug/ft3) 

0.002378 

0.002378 

0.001496 

0.000736 

0.002378 

0.000736 

Vm  (ft/sac) 

279.3 

1004.6 

1164.5 

584 

171.9 

2429.2 

q  *  PVa2/2(lb/ft2) 

92.7o 

1200 

1010 

126 

35-2 

2170 

ao  (dag) 

11.2 

2.1 

2.9 

7.5 

7-5 

0.8 

0o  (ft/sac) 

274 

1004 

1163 

579 

170.5 

2429 

WQ  -(ft/sec) 

54.2 

36.0 

58.9 

76.2 

22.4 

33.9 

(deg) 

0 

-7-4 

-3.8 

-*+•95 

-5.5 

-5.5 

-1.8 

-0.885537 

-0.779997 

-0.905967 

-0.714662 

-0.900031 

-0.643002 

-0.014978 

-0.013005 

-0.004499 

-0.022528 

-0.006670 

-0.004640 

C.235050 

0.094996 

0.039988 

0.225283 

-0.282005 

0.062601 

\ 

-0.568179 

-0.480165 

-0.393634 

-0.624148 

-0.600101 

-0.150010 

CBfta 

-0.855987 

-0.710445 

-0.63lt982 

-0.950073 

-1.14000 

-1.31291 

-0.660375 

-0.839857 

1.197740 

-0.820280 

-1.00135 

2.00191 

^q 

-3.787501 

-3-542457 

-4.056669 

-4.140056 

-6.12003 

-4.99989 

n 

-0.123596 

-0.079564 

-0.070374 

-0.113761 

-0.164400 

-0.006850 

Cr 

0.026750 

0.012195 

0.007703 

0.027221 

0.007410 

0.016651 

L6r 

% 

-0.289997 

-0.393388 

-0.405630 

-0.312417 

-0.441001 

-0.181412 

tt6. 

0.053330 

0.019453 

0.011381 

0.058537 

-0.077100 

0.009411 

ar 

0.184332 

0.069535 

0.063858 

0.147335 

0.500020 

0.032899 

CNa 

0.090368 

0.083167 

0.114795 

0.075179 

0.126700 

0.021700 

c*6r 

-0.101924 

-0.045250 

-0.023636 

-0.099188 

-0.123200 

-0.015897 

C"b. 

0.008559 

0.001806 

0.002284 

0.011366 

-0.017000 

0.001655 

c"p 

-0.000126 

0.020995 

0.040992 

-0.021406 

-0.012700 

-0.002981 

CKr 

-0.304765 

-0.322117 

-0.377810 

-0.305616 

-0.668000 

-0.146700 

a  *  375  ft  b  •-  38.7  ft  o  =  10.8  ft 

V  *  21,889  lb  M  =  680  slugs  eg  at  30  parcant  M3C 

Ix  *  13,635  slug-ft2  I  =  58,966  slug-ft2  X*  =  67,560  slug-ft2  IM  =  2,933  slug-ft2 
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